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Abstract. A generalized Cho-Faddeev-Niemi ansatz for SU(3) Yang-Mills is investigated. The correspond- 
ing classical field equations are solved for its simplest parametrization. From these solutions it is possible to 
define a confining central non-relativistic potential used to study heavy quarkonia. The associated spectra 
reproduces the experimental spectra with an error of less than 3% for charmonium and 1% for bottomo- 
nium. Moreover, the recently discovered new charmonium states can be accomodate in the spectra, keeping 
the same level of precision. The leptonic width show good agreement with the recent measurements. The 
charmonium and bottomonium El electromagnetic transitions widths are computed and compared with 
the experimental values. 
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1 Introduction 

Charmonium and bottomonium are essentially non rel- 
ativistic systems. Their dynamics can be understood in 
terms of a potential. Ideally, we would be able to derive 
such a potential directly from QCD. However, the con- 
fining potentials currently used to describe such systems, 
although being related to the fundamental theory, are not 
directly derived from QCD. 

In this paper, we report on a class of classical solutions 
of the Yang-Mills theory, computed after the introduction 
of a generalized Cho-Faddeev-Niemi |2*||^] ansatz for the 
gluon field. The solutions allow for a definition of a confin- 
ing potential which is then used to study charmonium and 
bottomonium spectra, their leptonic decays and charmo- 
nium electromagnetic transitions. For interquark distance 
between 0.2 fm to 1 fm, the new potential is essentially the 
lattice singlet potential. For larger distances it grows expo- 
nentially with the interquark distance and for smaller dis- 
tances is coulomb like. In what concerns the spectra, the 
new potential is able to reproduce charmonium and bot- 
tomonium spectra with an error of less than 3% and 1%, 
respectively, including the recently discovered new char- 
monium states X(3872), X(3940), F(3940), Z(3940) and 
Y(4260). If the spectra is reasonable well described, the 
hyperfine splittings 3 Si — 1 So turn out to be to small. This 
can be due to a missing scalar confining potential or the 
need for including contributions from other configurations, 
namely those related with the short distance behaviour of 
the theory. In what concerns the leptonic decays and the 
charmonium electromagnetic transitions, the new poten- 
tial is able to describe well the experimental data. Part of 
this work is described in pQ. 



2 Classical Configuration and a Heavy Quark 
Potential 

In the following we use the notation of 0]. Following a 
procedure suggest in , we introduce a real covariant con- 
stant field n a , 



D^n a = d^n a + gfabcn b A c u = . 



(1) 



Multiplying this equation by n a it cames out that n a n a 
can be chosen to be one. Defining the color projected field 



C M such that 



Al = n a C^ + X" n a X a = 0, 



(2) 



and replacing this definition in QJ, Xjl can De related to 
n a - see rjQ for details. For the simplest parametrization of 
n a = 5 al sin9 + S a2 cos 9 and for the simplest gluon con- 
figuration, the classical field equations in Landau gauge 
becomes 



d^e = o, 



(3) 



with gA^ = -5 a3 d^9. Note that there are no boundary 
conditions for 9. Equation J2J can be solved by the usual 
method of separation of variables. The solution considered 
here being 



A 3 = A(e At -b T e- At )v (r), 



A 3 = -(e M + b T e- M )vV (r), 
V (r) = A- 



( sinh(ylr) 



-AR 



(4) 
(5) 

(G) 
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Table 1. Charmonium Spectra. The right columns report the 
difference, in MeV, between theory and experimental numbers. 



Table 2. Bottomonium spectra. The right columns report the 
difference, in MeV, between theory and experimental numbers. 



Theo 
(MeV) 


Exp 
(MeV) 




Theo 
(MeV) 


Exp 
(MeV) 




J PC = 

3075 

3632 

4135 

4638 


o-+ 

Vc (lS) 
Y(4260) 


95 
-5 
-125 


[IP] 

3556 
3462 
3372 
3478 


Xc2(LP) 

Xei(lP) 
Xco(lP) 
h c {lP) 



-49 
-43 
-48 


J PC = 
3097 
3659 
4164 


I" [ 3 S 1 ] 
J/rl>(lS) 
i,(2S) 




-27 


J PC = 

3688 
4155 


I" [ 3 Pi] 
^(3770) 


-83 



M[J/iKlS)-r h (lS)] = 22\ Tbso : 
M[V>(2S)-^(2S)] =26| Theo 



117| Exp MeV 
48| Ex p MeV 



Assuming that the spatial part of Aq can be identified with 
a non relativistic potential for heavy quarkonia (see for 
details), the squared difference between V and the lattice 
singlet potential integrated between 0.2-1 fm was mini- 
mized to defined the various parameters; for r in MeV -1 , 
A = 11.254, B = -0.701, A = 228.026 MeV. The difference 
between Vb and the singlet potential is less than 50 MeV 
for r 6 [0.2 — 1] fm. In the following, the spin-dependent 
potential assumes that the interaction is pure vectorial 
type and that Vq describes the interaction between the 
heavy quarks. 



2.1 Charmonium and Bottomonium Spectra 

The charm quark mass, m c = 1870 MeV, was adjust to re- 
produce the experimental value of the Xc2(^P) — J/ijj(lS) 
mass difference. The bottom quark mass, m& = 4185.25 
MeV, was defined to reproduce the experimental value for 
the T(2S) — T(1S) mass difference. 

In table [I] we report the low lying charmonium states. 
Note that the masses were shifted to reproduce the J/ip(lS) 
experimental value. The third columns of table gives, in 
MeV, the differences between the theoretical prediction 
and the experimental values taken from 0. In what con- 
cerns the spectra, there is good agreement between the 
theoretical spectra and the measured values. The only 
particles which don't fit well in spectra are 0(4040) and 
-0(4415). In what concerns these two states, the experi- 
mental information is quite scarse and the particle data 
book comments that the "interpretation of these states as 
a single resonance is unclear because of the expectation of 
substantial threshold effects in this energy region" . Curi- 
ously, both particle masses are essentially the sum of J /ip 
with J PC = ++ known mesons. 

In what concerns the new charmonium states, the the- 
oretical predicitions which are close to the experimental 
values are: Z(3930), mass of 3929 ± 5 ± 2 MeV, J = 2, is a 
1 3 F 2 state with 3932 MeV or a 2 3 P 2 state with 4048 MeV 
; AT (3940), mass of 3942 ±11 ±13 MeV, and F(3940), mass 
of 3943 ± 11 ± 13 MeV, can be any of the following states 
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Exp 




Theo 


Exp 




(MeV) 


(MeV) 




(MeV) 


(MeV) 




J PC =0 


-+ 




J PC = 


I" [ 3 Si] 




9457 


Vb(lS) 


157 


9460 


r(is) 





10018 






10023 


r(2S) 





10380 






10385 


r(3S) 


30 


10721 






10727 


T(4S) 


148 


11059 






11065 


r(iio20) 


46 


[IP] 






[2P] 






9983 


Xm(1P) 


71 


10326 


X62(2P) 


57 


9941 


Xw(lP) 


48 


10283 


Xm(2P) 


28 


9894 


XW(1P) 


35 


10234 


Xw(2P) 


2 


9955 


h b (lP) 




10296 


h b (2P) 





[ 3 £l] 

10159 



10796 T(10860) 



6!) 



10476 r(4S) (?) 103 



10179 T(ID) 
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2 3 P , 3831 MeV, 2 3 Pi, 3938 MeV, or 2 1 P 1 , 3959 MeV; 
V(4260), mass of 4260 MeV, can be any of the follow- 
ing states 3 3 5!, 4164, 2 1 P> 2 , 4230 MeV, 2 3 P> 2 , 4230 MeV, 
1 1 G 4 , 4253 MeV, or 1 3 G 4 , 4260 MeV. For charmonium 
X(3872), mass of 3871. 2 ±0.5 MeV, the experimental data 
favors the following quantum numbers J PC = l ++ ,2 
For the potential considered here, the closest states with 
such quantum numbers are 2 3 Pi, 3938 MeV, and 2 3 Po, 
3832 MeV. 

In table|3we report the low lying bottomonium states. 
The masses were shifted to reproduce the Y(IS) experi- 
mental value. 

In conclusion, the non-relativistic potential Vq is able 
to explain the charmonium spectra with an error of less 
than 3% (~ 100 MeV) and the bottomonium spectra with 
an error of less than 1.5% (~ 100 MeV). 



2.2 Leptonic Decays of Charmonium and 
Bottomonium 

For the leptonic decays we follow the van Royen-Weisskopf 
approach and assume that QCD corrections factorize 
in the calculation of the widths, i.e. 



( 3 ^l) 



~mT- 1^(0)1 

25 e\ a 2 
2m\ M 2 



4 2 i>(o) 



(7) 
(8) 



times QCD corrections; e q is the quark electric charge, a 
the fine structure constant, m q the quark mass, M the 
meson mass, R n s(0) the S-wave meson radial wave func- 
tion at the origin and i^M ) the D 

-wave meson second 
derivative of the radial wave function at the origin. To 
avoid the problem of estimating the QCD corrections, in 
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Table 3. Charmonium and Bottomonium Leptonic Widths 
in KeV. The experimental figures are from [Sj- The limit for 
Y(4260) is from 0. 







Theo 


Exp 


Theo 


Exp 


2 3 Si 


i>{2S) 


2.84 


2.48(6) 


T(2S) 0.426 


0.612(11) 


3 3 Si 


Y(4260) 


2.16 


< 0.40 


T(3S) 0.356 


0.443(8) 


4 3 Si 








T(45) 0.335 


0.272(29) 


5 3 Si 








Y(5S) 0.311 


0.130(30) 



Table 4. El electromagnetic charmonium widths in KeV. In 
the calculation, for the meson mass it was used the experimen- 
tal measured mass. Experimental numbers are from 







Theo 


Exp 


i>(2S) — . 


Xc2(lP)+7 


30 


27 ±2 




Xel(lP)+7 


43 


29 ±2 




Xco(lP)+7 


51 


31 ±2 


Xc 2 (lP) — 


J/V(lS)+7 


414 


416 ± 32 


Xd(LP) — » 


J/^(15)+ 7 


308 


317 ±25 


Xco(lP) — » 


J/V(LS)+7 


146 


135 ± 15 


V'(3770)(l 3 Di) - 


- Xc2(lP)+7 


4 


< 21 




Xd(lP)+7 


110 


75 ±18 




Xco(lP)+7 


359 


172 ± 30 



tablc|3we report the theoretical witdhs computed relative 
to J/ip(XS) for charmonium and T(1S) for bottomonium. 

For charmonium, r e + e - for ip(2S) is slightly larger 
than the experimental value. For bottomonium, the witdh 
for the T(2S) is below the experimental number but the 
witdhs for Y(3S) and T(4S) agree with the experimen- 
tal figure within two standard deviations. For T(5S), the 
theoretical prediction is a factor of 3 higher than the ex- 
perimental figure. For D-wave mesons, the theoretical pre- 
dictions for the leptonic widths are much smaller than the 
measured _T e + e -. However, for cc if one considers mixing 
betweem S and D mesons, it is possible to bring all widths 
quite close to the experimental figures; see £Q for details. 



Table 5. El electromagnetic charmonium widths, in KeV, for 
the new charmonium states. In the calculation, for the me- 
son mass it was used the experimental measured mass. Only 
transitions to wellknown states are considered. 









r 


X (3872) [2 3 Pi] — > 


JM1S)± 7 
^(2S)± 7 


697 
181 


48 
63 


A" (3940), Y(3940) 
2 3 Pi or 2 3 P — ► 

2 3 Po — ► 
2 x Pi — ► 


JMis) + 7 

V>(2S)± 7 
^(3770) ±7 
Vc(lS) ±7 


755 
249 
168 
845 
293 


61 
165 
70 
86 
271 


Z(3930)[2 3 P 2 ] — > 


J/V(1S)± 7 
V(2S)± 7 


744 
235 


59 
140 



Table 6. El electromagnetic bottomonium widths in KeV. 
The experimental numbers are from [10| . 







Theo 


Exp 


T(2S) — ► 


Xc 2 (lP)±7 


3.30 ± 0.04 


2.21 ±0.16 




Xci(lP)±7 


3.21 ±0.04 


2.11 ± 0.16 




Xco(lP)±7 


2.10 ±0.02 


1.14 ±0.16 


r(3S) — > 


Xc2(2P)± 7 


3.22 


2.95 ±0.21 




Xci(2P)± 7 


2.95 


2.71 ± 0.20 




Xc0 (2P)±7 


1.82 


1.26 ± 0.14 



To finish, in table we report the bottomonium El 
electromagnetic transitions. 
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2.3 El Electromagnetic Transitions 

For the computation of the El electromagnetic charmo- 
nium transitions we follow |S|. Tabled report the theoret- 
ical estimates of the decay widths known experimentaly. 
On overall, the agreement between theory and experiment 
is good. The exception being the transitions involving the 
scalar meson x c o(lP). Remember that for the meson spec- 
tra, the larger deviations from the experimental numbers 
occured for the scalar mesons. 

Given the good agreement for the known experimental 
radiative transitions of the charmonium, for the theoreti- 
cal states which can describe the new charmonium states 
we have computed their El electromagnetic widths. In ta- 
ble we report those widths which are larger than ~ 40 
KeV. In principle, the electromagnetic transitions can al- 
low to distinguish not only the various theoretical models 
but also the difference meson states. 
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